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SUMMARY 

The development of rotary wing induced-velocity theory i? traced from its 
origin as a momentum-theory estimate of average interference, vhrough simple 
vortex theory, to its present status where it is indispensible in calculating 
blade loads. Applications to a v\riety of interference problems is demonstra- 
ted. Wherever possible, experimental results are presented to confirm the 
theory. 


INIEODUCTmN 

All aerodynamic lifting devices, including helicopter and autogyro rotors, 
obtain their lift as the reaction from pushing air downward. The momentum thus 
imparted to the air sets up an induced field which interacts or interferes to a 
significant degree with anything placed near the lifting system; indeed, it 
interferes with the lifting system itself. 

In the case of a helicopter rotor, the velocities induced by the field 
influence the forces on the blades as they rotate; they alter the effectiveness 
of tail surfaces; and they have powerful influences on the interference with 
the additional lifting devices used on a multitude of converti plane designs. 

The detail in which the induced flow must be known veries with the intended 
application. A gross overall average velocity through the rotor is adequate for 
many performance problems. The special, time-averaged, field is often adequate 
to predict tail effectiveness or mutual interference with other lifting systems. 
If detailed load distributions over the blades are required, an equally detailed 
knowledge of the instantaneous induced velocities along the blades is required. 
The difficulty encountered in obtaining the required infonnation is more than 
proporticn.al to the detail required in the velocities. In the simplest cases, 
the calculation is in closed Form and occupies three or four lines. In the more 
difficult cases, the calculations may exceed the capacity of some of the largest 
digital computers. 

The present paper attempts to summarize the state of knowledge of rotary- 
wing flow fields. It begins with the simple momentum methods applied to the 
rotor by tlauert. It then turns to the time-averaged flow field, where the 
first truly usable results were obtained by Walter Castles of Georgia Tech. 

The extensicn of Castle's work by NACA and a comparison of the theory with flow 
measurements are presented. Finally, the modern efforts toward using more 
detailed d’’qital methods to obtain blade-load distributions are described. 
















In most cases, figures and equations presented herein, have been abstracted 
from the original reports. The notation, which varied from author to author, is, 
therefore, not consistent throughout this paper. It is hoped that the notes 
presented with equations, and the general context of the discussion, will suffice 
to allow the reader to infer the proper meanings for the symbols throughout the 
development. 


DISCUSSION 

MOMENTUM THEORY - THE AVERAGE INDUCED VELOCITY 

Glauert (ref. 1), in the early 1920's, began the development of rotary wing 
blade-element theory. It was clear from analogous work on propellers and wings 
that the local angles of attack could not be found solely from geometric con- 
siderations. There must be an interference, or induced, velocity which resulted 
from the lift, and which entered into the angle of attack equations. No theory 
of the time predicted this induced velocity; however, the Fro:d momentum theory 
for propellers in static thrust resulted in 

T = -pttR Wq (2wq) 

where T is the thrust, p the mass density of the fluid, R the radius, and 
Wf the induced velocity through the disk (positive upward for a rotor). On the 
other hand, for a wing having a span of 2R, it was known that the lift L was 
given by 

L * -pttR^V(2Wq) (2) 

Glauert noted that these two expressions could be combined in a continuous 
fashion if the total velocity through the rotor disk v/as taken as the vector 
sum of the induced and forward velocities; that is, if T ^ L 

T = -p7rR^Vp(2wQ) (3) 

where is the absolute value of the sum of the induced and forward velocities. 
He applied the resulting value of wg as a uniform induced velocity, or down- 
wash, over the entire rotor disk. In terms of present rotor notation, where the 
thrust coefficient is 

C = T 

^ PttR^ (QR)^ (^) 


the inflow ratio (velocity ratio normal to the disk) is 


X = 


Vsina + Wq 
HR 


(5) 
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and the advance ratio (velocity ratio parallel to the disk) is 


_ Vcos a 


( 6 ) 


equation (3) may be expressed as 
, C^QR 

““ ° 

Here the matter rested for about 25 years. Within the usable performance 
range of the machines of those days, such a simple theory was reasonably ade- 
quate. The first advances in rotor flow-field theory were not really occasioned 
by great known inadequacies in rotor theory as such - rather, it was the need to 
estimate downwash angles at the wings and tails of unloaded rotor configurations. 
For these purposes, simple momentum concepts were totally inadequate, and it was 
necessary to turn to at least a rudimentary vortex theory or its equivalent. 


SIMPLE VORTEX THEORY 

Work on a simple vortex theory for rotors was begun during World War II by 
Coleman, Feingold, and Stempin (ref. 2); however, it was an additional ten years 
before this work was brought into ful'' flower by Castles and De Leeuw (ref. 3). 
(A mathematically different, but physically identical study was done in England 
by Mangier and Squire (ref. 4)). 


The Vortex Model 

The first step in setting up a simple vortex model is to choose a vortex 
pattern that is reasonably adequate on a physical basis and, at the same time, 
sufficiently simple that numerical results can be extracted. The wake of a 
rotor is more complicated than that of a wing, and a more involved model is 
required. If one examines a single-bladed rotor, as in figure 1(a), there 
would be a bound vortex along the blade. At the root, the free vortex is 
carried downward by the flow through the rotor and rearward by the flow past 
the rotor. The tip vortex is carried off in the same manner, but it also is 
affected by the rotation of the rotor; thus, it lies on a spiral path as shown. 

Even this simple spiral Is difficult to deal with so the next step is to 
divide the wake vortices into their axial and circumferential Cu„iponents and to 
discard the axial components as being of secondary importance. The wake then 
appears as a skewed stack of vortex rings as in figure 1(b). The final step is 
to assume that the rings are spaced so closely that they ai'e equivalent to a 
continuous elliptic cylinder of vorticity as shown in figure 1(c). 
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The Induced Velocities 


The induced velocities constituting the »*otor field may then be found by 
integrating the Biot-Savart law over the t tire wake. This can be acconplished 
in several manners. The development presented herein follows that of reference 
5. The induced velocity of an element in the wake (fig. 2) is 


dq = 


JL dr 

4it 


dsXa 

lal^ 


dL 


(8) 


where q is the induced velocity vector, L is the length along the wake, and 
Dr/dL is the vortex density (or vorticity) along the edge of the wake. 

From figure 2, s and a, are found to be 

i = T(R.'Os'^'+LsinX) + j(Rsin4') + M-LcosX) 


a = 1 (Rcos'^ + LsinX-x) + j(Rsin»P-y) + k(-Lcosx-z) 


and ds is found by differentiating s with respect to 
ds = [i(-sin>i') + j(cos>i^) + R(0)jRd'i' 

Substituting these values into equation (8) yields 


q = 


dr 

Jl 


2tt 


T j R 

-sin*^ cos'^' 0 

(Rcos'i'-Hsi nX-z ) (Rs i n'^-y ) ( - Lcosx-z ) 


0 0 


'(Rcos'^’+LsinX-z)^ + (Rsini/'-y)^ + (-LcosX~z)^ 


dLdii' 

T~- “ 


(9) 


w = 


For the moment, we will be interested only in the k (or w) component of 
induced velocity since this component is perpendicular to the disk. From 
equation (9) 


4it 



(xcos'i'+ysinij^-R-LsinxcosiiO dLdi|; 




R^+x^+y^+z^-2R(xcos'l'+ysin’/’) + 2L(zcosX-xsinX+RsinXcos’J0 + L*" 
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roj col 



The integration with respect to L is not difficult and the forms will be 
found in tables of integrals. After substituting limits and the nondimensional- 
izing, the result is 


w = 


-1 dr 
W “3l 


2tt 



i cos*. ^ 


sin'i' 


^ T 


sinxcos^' 


-FT + COS'/' 


L R 


sinX+ p cosX 


d'/' 


( 11 ) 


where R , the distance from the point (x,y,z) to the edge of the disk at 
is 


“ V 


R^+x^+y^+z^-2R(xcos'i'+ysin'/') 


(12) 


Equation (11) is not integrable, in general, except in terms of elliptic 
integrals of the third kind, and it is simpler to evaluate it numerically in a 
digital computer. One case which can be integrated is wh-;n (x,y,z) is at the 
center of the rotor (y = y = z = 0, R(. = *^) for which equation (11) reduces to 

2tt 

w = 1 dr 

'^0 ~3i- 

0 

The nondimensional induced velocity is obtained by dividing equation 
(11) by equation (12) to yield 


di/) = - 


1 dr 
7 ar 



2i\ 



R ' 

1 - cos'/'+ sin'l'j+ sinXcos’/'j 


fP • 

Lt* (cos*- si ox. I COSX 


R 

c 


dip 


(13) 


Local values of the induced velocity ratios may be obtained by numerical 
evaluation of equation (13). Charts of the induced velocities are available 
from which approximate values may be obtained. A collection of such charts is 
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given in reference 6. A sample set of charts (from ref. 7) for both the 
vertical and longitudinal induced velocity ratios is given in figure 3. 


The Average Induced Velocity 

As noted earlier, Glauert's formulation of the induced ’'elocity of a rotor 
was merely a plausible guess. The simple vo<ex theory can be used to examine 
the validity of the value which he selected. Consider a rotor of b blades, 
each of which has a bound circulation r over its entire length. The 
elemental lift on the blade is found from the Kutta-Joukowski theorem to be 

dT = pU.^rdr (^4) 

where Uj is the tangential velocity normal to the blade axis and r is the 

bound circulation. The total thrust of the rotor' is found by integrating the 
elemental thrust over the length (radius) of the blades and then averaging this 
value around the disk; that is 

2it R 

p (P.riMi^Rsin'l’) I'drd'l' = ^ bf^'^R^r (15^ 



Substitute equation (15) into equation (4) to obtain 
C = br 
''' 27tR^P 

and then solve equation (16) for T to yield 
2nirR^CT 

r = — ^ — L_ 
b 


(16) 


(17) 


Now the number of "vortex rings" formed per unit time at the rotor is 
BP/2 tt . These are carried down the wake with the mean velocity at the rotor 

(W + X^) . Thus, the vorticity along the edge of the wake is 


dr 

au 



/ 2 2 
PR y + X^ 


(18) 


Substitute equation (17) into equation (18), and simplify the result to >ield 


dr 

ar 


PRCt 


"2 ? 

' ' + X 


(19) 
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Simple vortex theory gives the induced velocity at the center of the rotor 
as equation (12). Substitute equation (19) into equation (12) to obtain 





( 20 ) 


Equation (20) for the center of the rotor is identical to Glauert's result 
(eq (7)) ^or the average over the entire rotor. The induced velocity varies 
widely o.^r the rotor disk; however, reference 8 has shown that, even when 
harmonic variations in r occur at the blades, the value at the center is 
identically equal to the average value. Thus, the simple vortex theory confirms 
the momentum theory formulated by Glauert. 


VORTEX THEORY FOR 

NONUNIFORM AXISYMMETRIC DISK LOADING 

At this point, only blades of constant circulation r have been considered. 
For such a blade, the radial gradient of blade loading is obtained by taking 
Uy = fir and reforming equation (14) to obtain 

§ = pfirr (21) 

Equation (21) shows that, for unifonn bound circulation, the thrust per unit 
length of blade varies linearly from zero at the root to a maximum at the tip. 

Now consider the thrust per unit area of the swept disk. The area of an annulus 
of radial width dr is 27 t rdr. Thus, the load per unit area on this annulus 
is (for b blades) 

dT _ bpfirrdr _ bpfiF /oo\ 

dA 2iTrdr 2 tt 

Thus, uniform circulation along the blades produces a uniform load per unit 
area on the swept disk. Because, the disk is circular, the spanwise distribucion 
of lift viewed from behind the rotor is semicircular; that is, a special case of 
an elliptic spanwise load distribution - the ideal distribution for the minimum 
induced drag of a wing. 

In general, the disk load distribution of a rotor is not uniform. If the 
radial distribution of induced velocity was neglected, a constant-chord untwisted 
blade would have a blade loading proportional to r^ not r as in equation 20. 
Twist and taper further alter the load grading, and the cutout to provide the 
rotor hub is a nonuniformi ty in itself. Any practical rotor has zero disk-load 
distribution at its center. 

In general, r is some function of r. Regardless of the actual variation 
chosen, equations (21) and (22) show that the disk load distribution will vary as 
1/r times the blade load distribution. The flow field of nonuniformly loaded 
rotor disks can be obtained from the field of the uniformly loaded rotor disk 
provided that the disk load distribution required is axisymmetric (ref. 9). 


Consider the triangular load distribution of figure 4. It can be 
approximated by a stepwise distribution of loads comprised of one large positive 
load (1.5 times normal strength) and a series of small negative loads of smaller 
radius. Such a load distribution would be generated by a nest of coaxial vortex 
cylinders such as shown in figure 5. 

The contribution of each of these cylinders to the flow at (X/R, Y/R, Z/R) 
will be proportional to its vortex strength and can be read from charts such as 
figure 3 at the ooint [ (X/R)(R/Rv) , (Y/R)(R/Rv), (Z/R)(R/Rv)] • The induced 
velocity at that point is the sum of the contributions of all the cylinders. 

A sample calculation for a triangular loading at (X/R = 0.2, Y/R = 0, 

Z/R =0.1) and for x * 90° is shown in figure 4. 

The actual distribution of load on the disk has very large and significant 
effects on the distribution of induced velocity. A few sample charts are pre- 
sented in figure 6. A direct comparison is presented in figure 7. Obse’^ve 
that there is a region of upwash in front of, and encompassing, the leading 
edge of the disk in each case. The major differences occur over the disk and 
across the wake. In constrast to the smooth increase in downwash along the 
rotor axis in the uniform case, the triangular case shows an increase at first, 
then a decrease to zero at the rotor center, an upwash behind the center, and 
then, finally a rapid increase to the trailing edge. The distribution across 
the wake itself is markedly different in bnth cases. 

One other point is very significant. Under the present assumptions, the 
rotor-induced flow is dependent only upon the skew angle x. The total flow 
field is the vectorial sum of the uniform external flow and the rotor induced 
velocities. The total field may be markedly different dependent upon the rotor 
angle of attack and loading. This feature is illustrated by figure 8. 


COMPARISON WITH EXPERIMENT 

NACA Rep 1319 (ref. 9) presents the results of comprehensive flow surveys 
near a rotor operating in the wind tunnel. Sample results are given in figures 
9 and 10. Comparison of theory and experiment in the central plane of the rotor 
is shown in figure 11. It is evident that the radial distribution of loading is 
significant in the correlation. On the other hand, an earlier paper (ref. 10) 
showed that the uniforri load distribution does give a reasonably good approxi- 
mation to the average induced velocity across a span equal to the rotor diameter. 

Another significant feature of figure 11 is the evident failure of the 
theory in the region aft of the rotor. This resjlt is obtained because of the 
extremely rapid roll up of the wake vorticity behind the rotor. Figure 12 shows 
contours of equal vorticity behind the rotor and it is evident that the roll up 
is already well underway at the trailing edge of the rotor. In this region, the 
wake can be modeled by assuming it to be about the same as that behind a wing of 
equal span. 
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Some Applications of the Theory 

Wind tunnel flow measurements give reas'^n to believe that the major ."actors 
affecting the flow field are included in the theoretical considerations desc'^^bed 
above. Thus, the theoretical results should be adequate to explain the induced 
effects found on helicopters and compound aircraft (ref. 11). This section of 
this paper will examine the extent to which this conclusion is true. 

Multi rotor helicopters .- The first example will be a tandem-rotor system 
which was tested in the Langley full-scale tunnel (ref. 12). Basically the 
model was simply two identical 15-foot-diameter rotors with no vertical offset 
and a gap of only 6 inches between the rotor tips. Thrust and torque were 
measured and, in addition, flow surveys (ref. 10) were made at locations very 
near to the rotor. For one level-flight condition, the flow measurements are 
shown in figure 13. The data points represent th', average induced velocity 
across a span of one rotor diameter. These points indicate that, the average 
induced velocity through the rear rotor was three times that through the front 
rotor. This correspond? with the simultaneous power measurements which indi- 
cated that the rear rotor required three times as much induced powe- a" the 
front rotor. 

The calculated induced velocities in this case are also shown on figure 13. 
These calculated curves are very close to the measured flow data and indicate 
that the rear rotor should require more than 2.8 times as much induced power as 
the front rotor. Thus it appears that the theory should be adequate for esti- 
mating the effect of interference on the overall performance of multi rotor 
helicopters. There are, however, certain reservations to this conclusion which 
will be discussed in a later section. 

Rot or-wing interference .- Recent attempts to combine the forward-flight 
capablTTty of an airplane with the vertical -flight capability of a helicopter 
have led to a number of un'ioaded-rotor convert! plane designs. Generally 
speaking, these compound aircraft are characterized by ine presence of a rotor 
and a wing operating in close proximity to e :h other. Some degree of mutual 
interference is bound to be present and it may have a major effect upon the 
overall performance of the aircraft. 

Two such machines are shown in figure 14. These are the PCA-2 autogyro and 
the XV-1 convertiplane. If the hovering ability of the XV-1 Is neglected, the 
two aircraft are remarkably similar despite the passage of some twenty-five years 
between their designs. No power is applied to the rotor in forward flight since 
both operate as autogyros in this flight regime. An auxiliary wing is provided, 
and this wing carries a large portion of the lift. The propulsive force is 
obtained in each case by means of a propeller. The major differences between 
the two aircraft lie in both the absolute and the relative loadings of the wing 
and the rotor. 

The PCA-2 autogyro was tested in flight. The wing load was determined by 
measuring the pressures on the wing surfaces. The remainder of the gross weight 
is supported by the rotor. Thus the performance is relatively simple to compute 
since the rotor lift coefficient and, consequently, its mean, or momentum theory. 
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volue of Induced velocity are already known. The induced velocity at the wing 
due to the rotor is then taker, from the charts. Adding this induced velocity 
vectorially to the free-stream velocity yields the direction of flow at the 
wing and the effective wing angle of attack, ‘’’he lift of the wing can then be 
detertnined from the isolated wing characteristics. 

Figure 15 compares the wing load in percent of the total lift as 
calculated by the foregoing procedure and as measured in gliding flight. (If 
the powered-flight conditions had been chosen it would, of course, have been 
necessary to correct the wing lift for the increased dynamic pressure in the 
slipstream propeller.) At each tip-speed ratio, the measured load is repre- 
sented by a circle and the calculated load is represented by a square. The 
agreement is good, the maximum error being about 6 percent. 

The XV-1 converti plane was tested in the Ames full-scale tunnel and, 
because of the manner in which the tests were conducted, presents a much more 
interesting problem. During the tests, the performance of the rotor, the con- 
verti plane without its rotor, and the complete converti plane were all measured 
separately. The results are shown in figure 16, where the symbols indicate the 
measured ^ift and drag coefficients which are based in all cases on the rotor 
disk area. The sum of the "rotor only" data and the "rotor off" data is given 
the broken line. The difference between this line and the data for the 
complete configuration is large and may be charged almost entirely to 
interference. 

The calculation in this case depends upon a knowledge of the total lift 
required, just as it does in a flight calculation. At each angle of attack, the 
rotor is assumed to have its "isolated" lift. Then the wake skew angle is com- 
puted and the induced velocity at the wing is determined. This induced velocity 
is used to determine the effective angle o^ attack and lift of the wing. The 
sum of rotor lift and wing lift, in general, is not the correct lift. The lift 
of the rotor is altered to correct the overall configuration lift, and the pro- 
cedure is repeated iteratively until the total configuration lift is obtained. 

The process is rapidly convergent in a practical sense. Once the lift coeffi- 
cients of both the rotor and the wing are determined, it is a simple matter to 
determine thei*" drag coefficients. 

The results of such calculations are shewn as the solid curves of figure 16. 
The accuracy of the procedure is attested by the remarkably close agreement with 
the observed perfcrmance of the complete configuration. 

Tail Surfaces .- Figure 17 presents the downwash measured from the floating 
tail of the XV-1 . Similar measurements obtained from an 8/31 scale model in 
the return passage of the WADC 20-foot tunnel are also shown. 

The induced angle at the tail is affected by both the wing and the rotor. 

The lift coefficients obtained from the rotor-wing interference calculations 
were used for determining the contribution of each to the downwash angle at the 
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tail. The calculated results are shown in figure 17 for the cases of uniform 
and triangulat' loading. It should be noted that the span of t^e tail is approx- 
imately one-half radius. Thus this case corresponds to neithet the one-diameter 
span required by the uniformly loaded rotor induced-velocity chart nor the zero 
span required by the triangularly loaded rotor chart. Therefore, the data 
should fall between the two calculated curves, and tnis is exactly the result 
shown in figure 17. 

The degree of agreement shown here can be expected - ily when the rotor lift 
coefficient is low. If it were much higher the flow would roll up faster and 
might cause large errors. For the case of large lift coefficients, the downwash 
angles may be calculated fairly well be assuming that the flow is similar to 
that of the uniformly loaded wing. 

Aerodynamic loading of rotor blades .- One of the most troublesome features 
of rotary-wing aircraft is the large magnitude of the fluctuating aerodynamic 
loads on the rotor blades. These loads have been measured during wind-tunnel 
tests of rotors equipped with instrumentation capable of measuring the instan- 
taneous pressure distribution on the blades. Figure 18 shov/s the aerodynamic 
loading at the threa-quarter-radius station as measured during one of these 
tests (ref. 13). It may be seen that in this case the loading fluctuates with 
an amplitude of ubout one-half of the mean load. Needless to say, the stresses 
resulting from fluctuating loads of this magnitude could easily represent <-he 
critical design condition for the rotor. 

The data shown in figure 18 was obtained at a tip-speed ratio of only 0.08. 
At this speed, the older blade-element theory with uniform inflow indicates that 
the aerodynamic load should be essentially constant. This, in itself, is not 
too disturbing since it is generally recognized that the assumption of uniform 
induced velocity inherent in the theory limits i' 'o computing only the over? 
performance rather than the detailed blade loading. 

The blade-element theory should, of course, yield a much more accurate 
estimate of the blade loading if the correct induced velocity were used. In 
order to find just what this induced velocity should be, one investigation 
took the measured blade loading, inserted it into blade-element theory, and 
worked backwards through the theory to obtain an effective induced velocity. 

The results are shown in figure 19 where they are compared with the theoretical 
induced velocities. It is immediately evident that there is no correlation at 
all between the required induced velocity and the theoretical induced velocity. 
Even the trends are completely reversed in some locations. Thus it appears that 
the calculated induced velocities cannot lead to the correct blade loading. 

Reexamination of basic concepts .- The complete failure of the induced flow 
theory in the blade-load problem is in sharp contrast with its success in com- 
puting interference effects. The reason for this failure is very basic and may 
be explained by returning to the original concept of the rotor wake as a spiral 
vortex (fig. 20). As pointed out previously, this spiral will move downward 
with time and the actual induced flow will be a function of time. This time- 
dependency was lost as a direct result of the simplifications adopted for the 
analysis. The calculated induced velocity is therefore a time-averaged value 
of the actual induced velocity. 
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Now consider a point P , fixed in its position with respect to the rotor, 
as, for example, a point on a fixed-wing surface. As the spiral vortex pro- 
gresses downward, the point P will "see" all of its possible positions. 
Therefore, the mean induced velocity at P should be the time-avereged value 
calculated by theory. Thus, the calculated induced velocities allow accurate 
estimates of rotor-wing interference to be made. 

Next consider point Q which is on the rotor blade and which rotates with 
the blade. Since the spiral vortex starts at the blade tip, there will, at any 
azimuth position, be one, and only one, position of the spiral with respect to 
Q . Thus Q will not "see" the time-averaged induced velocity giv?n by the 
theory. Thus, as shown before, the theory cannot be applied to finding the 
induced velocities needed for calculating he blade loads. The actual blade 
loading problem is closely related to flutcer in that the unsteady aerodynamics 
cannot be neglected. 

The more modern rotor theories take into account the unsteady character of 
the flow and, in addition, turn to mo’/e refined models of the wake. A brief 
survey of some of the more recent developments in this area will folio'-;. 

There is an intermediate problem between the average external interference 
and the detailed blade loads. Consider the point R on a second rotor near the 
first. Does t.his point see a time-averaged velocity or not? Locally, on the 
blade, it probably does not, and the change in the detailed blade loading o'; a 
result of intertr.ence can probably not be calculated with any degree of 
accuracy. However, at least for the tandem rotors of NACA TN's 3236 (ref. 12) 
and 3242 (ref. 10), there is experimental evidence to indicate that the overall 
rotor forces can be determined by use of the present theory. This may not be 
the case for all configurations. It is therefore suggested that multirotor 
interference calculations based on the present induced-fl'^w theory be accepted 
for the present but always subject to the provision that experimental verifi- 
cation is required in each case. 


CALCULATION OF BLADE LOADS 

The 'ore modern applications of rotary wing vortex theory involve the 
calculation of the details of the blade load distribution. The discussion 
herein will largely follow the developments in the United States. Although 
complete details do not appear to be available in this country. Mil's textbook 
(ref. 14) makes it clear th:t very similar developments have taken place in the 
USSR and that those analyses often predated our ov;n. The present discussion 
must also deal only in generalities. Closed form solutions to the problem do 
not exist. Only the numerical results from specific sample cases run on digital 
computers are available. 

The earlier aiscussion has shown clearly that the instantaneous velocities, 
not the average velocities, at the blade are required. In addition, it is 
evident that quasi-static solutions can not be meaningful in the face of the 
widely varying local velocities and angles of attack on the blades. The need 
for an unsteady analysis was further shown by the essentially simultaneous papers 
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of Loewy (ref. 15) and Timmon and van de Vooren (ref. 16). Both analysis 
demonstrated the existence of previously uncalculated wake-excited flutter 
modes even for a hovering rotor. 

The initial theoretical investigation of blade loads was conducted by 
Piziali and DuWaldt at Cornell Aero Labs (ref. 17). The model used for their 
calculations is shown in figure 21. The blade is divided into a number of 
radial segments with the bound vortex fixed on the quarter chord. The bound 
circulation along each radial segment is constant. Thus, at each junction of 
segments a vortex will be shed, and the strength of this vortex will be equal 
to the difference in circulation of the two segments. The rotational motion 
of the blade is considered in terms of a number of azimuth steps with the bound 
circulation changir^ ^rom step to step. In order to satisfy conservation of 
vorticity, a vortex must be shed from, and parallel to, the trailing edge of the 
blade at each step. Furthermore, the strength of this vortex must be equal to 
the difference in bound circulation between the two azimuth steps at each radial 
station. This grid-like wake is assumed to be carried uniformly away from the 
rotor tip-path plane by the forward velocity and by a uniform average induced 
velocity. Thus, it is constrained to lie within the cylindrical boundaries of 
the wake model used earlier in the average downwash calculations. 

At the start of the calculations, the bound circulations on the blade are 
unknown; however, all of the wake vortex strengths may be formed as functions 
of the unknown bound circulations. The blade angles, flapping angles, and tne 
deformations of the blade were preassiqned in this analysis; thus, the dynamic 
response of the blades was not calculated. Control points were specified on 
the three-quarter chord at the midpoint of each segment. It was then possible 
to set up an array of simultaneous equations in the unknown blade circulations, 
using the Biot-Savert law to obtain the induced velocities at the control 
points. These simultaneous equations were then solved by matrix inversion tech- 
niques on the computer. 

Observe the interrelationships involved in tiie solution. Charts of induced 
velocities for standardized load distributions no longer have any meaning. The 
only distribution with meaning is the one for that particular rotor at that 
particular flight condition. 

The results from the n^w techniqu° v.'ere far superior to the old 
ur.iform-downwash, quas’-static theory in many respects; however, as pointed out 
by the authors, there were a, so numerous problems. One major problem was that 
the calculated first-harmonic flapping moment wa^ not zero, even when measured 
blade motions and deformationi; were input co the program. It was further 
recognized that the single three-quarter-chord control point was probably too 
crude, and that wake deformations might have to be considered as well. 

Certain of these improvements were incorporated by Piziali in a later paper 
(ref. 18). In this version of the analysis, the blade was represented as a 
continuous distribution of vorticity by means of a Glauert series anL the span- 
wise vorticity behind the trailing edge was "effectively fitted" with a con- 
tinuous sheet (ref. 19) in the region immediately behind the blade (fig. 22). 

The entire calculation was set up in a loop such that the airloads were computed 
and then the resulting flapping and blade bending were computed. These results 
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were compared with the same quantities from the previous iteration and the 
entire procedure was repeated until the airloads, flapping, and dynamic 
response all converged. A sample set of results is compared with flight 
measurements in figure 23. Considering the complexity of the calculation and 
the simplifying assumptions, the results were quite encouraging. 

Piziali's discussion pointed to the fact that many of the remaining 
discrepancies between theory and experiment could be qualitatively accounted 
for by physically reasonable distortions o^ the wake from the simple helical 
sheets which he was still using. Actually, Crimi , also at Cornell Aero Lab, was 
already working on the problem. 

Crimi 's original analysis (ref. 20) was directed at the external field of 
the rotor, and, for this purpose, he represented the wake by a single tip 
vortex from each blade No spanwise vorticity was used despite the assumption 
that the tip vorticies had a strength containing a sine component. Rigid 
helices were assumed at the start of the calculation, and numerous control 
points were chosen along the helices. Then the induced velocities at the 
control points were calculated. (Certain infinites in the calculation were 
avoided by the assumption of finite vortex cores.) The vortices were then 
allowed to convect to a new position with these velocities over a time period 
equal to one step in olade azimuth position. The process was repeated with the 
new positions until further changes in the wake positions were insignificant. 

The final wake position could then be used to calculate the flow throughout the 
field. 


The wake resulting from this calculation is shown in figure 24. It 
obviously duplicates, at least qualitatively, the wake distortions which were 
found earlier in the wind-tunnel flow measurements. A comparison of the time- 
averaged measured and calculated flow fields is shown in figure 25. Considering 
that Crimi 's assumptions omit the effect of radial disk-load distribution, the 
agreement is very encouraging. 

More recently, Sadler (refs. 20, 22) has incorporated similar wake 
distortions into the calculation of blade loads using a full vortex-mesh and 
working with the unsteady aerodynamic problem. The solution is "’ided by a 
"rotor start-up" approach; that is, the theoretical rotor commences from a dead 
start in which there is no wake at all. in the first time increment a short 
wake is sheu, approximate blade loads are calculated, and the short wake is 
relocated. The process is repeated by steps in aximuth until the complete wake 
is generated. Samples of the results are presented in figures 26 to 31. Once 
more the results are encouraging. 

It would be premature to claim that the problem of calculating blade loads 
is solved. The comparisons between theory and experiment still indicate 
uncomfortably large differences. On the other hand, it is evident that many of 
the necessary elements to a solution have been determined and the basic problems 
dealt with to some extent. It is clear that far more detiil is required than 
has been employed to date. Unfortunately, this detail is very expensive in terms 
of computer storage requirements, running time, and cost. 
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CONCLUDING REMARKS 


The development of rotary wing induced- velocity theory has been traced 
from its origin as a momentum- theory estimate of average interference, through 
simple vortex theory, to its present status where it is indispensible in 
calculating blade loads. Applications to a variety of interference problems 
has been demonstrated. Wherever possible, experimental results have been 
presented to confirm the theory. 
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Figure 3. - Induced velocities in the central longitudinal plane of 
rotor with uniform disk loading. 
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aure 3. - Continued 







Contours of u/w^ 

(d) X = 30° ° 
Figure 3. - Continued. 
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Figure 4. - Triangular disk 








(a) Uniform load ^ ' 

- ^ (b) Nonuniform axially syirmetric load. 

Figure 5. - Assumed vortex pattern of rotor wake 
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Figure 6. - Continued. 
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Figure 7. - Comparison of induced velocity ratios for uniform and 
triangular disk load distribution, x = 76o. 
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(a) Uniform disk loadi'ig. 

» 

Figure 8. - Flow field at different angles of attack h a single 

wake skew angle of 50° . 








Figure 9. - Measured stream angles near a rotor, y = 0.14 
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Figure 10. - Contours of dynamic pressure behind a rotor. 
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Figure 13. - Downwash of tandem rotors. 
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Figure 15. - Division of load between wing and rotor of PCA-2 autogyro. 



Figure 16. - Performance of XV- 1 at 100 knots. 





Downwash at tail 
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StOEMfHTEO TJUIUIIO VOUTICES 

SEaEMENTED SHED VOUTICES 

Figure 21. - Vortex wake used by Piziali and Duwaldt in reference 17. 



Figure 22. - Blade representation used in reference 18 by Daughaday and Piziali 
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Figure 25. - Comparison of computed^ and measured Induced 
downwash below the rotor plane, u = 0.14, C = 0.00371, 
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Figure 27. - Wake model with combination of "full mesh" wake 
and "modified" wake of trailing vortices only. 
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Figure 28. - Comparison of computed and measured induced 
below the rotor plane, u = 0.14, C = 0 00371. x = 0. 
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Figure 30. - Comparison of computed and measured induced 
downwash below the rotor plane, u = 0.14, C = 0.00371, 

X = 0.5. 









(a) plan view 



(c) rear view 


Figure 32. - Single two-bladed rotor's tip vortex locations 

y = 0.141 a. = -5.10. 
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Figure 33. - Blade airloads for H-34 helicopter, Ib/ft versus 

azimuth, y = 0.2 
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